Hole Binding around Ni Impurity in Cuprates 
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We examine the influence of Ni impurity in cuprates on the distribution of hole carriers by 
performing numerically exact diagonalization calculations for a model consisting of Cu3ci, Ni3ci, 
and 02p orbitals. Using realistic parameters for the system, we find that a hole is predominantly 
bound to 02p orbitals around the Ni impurity forming the Zhang-Rice doublet. This imposes strong 
restrictions on modeling Ni-substituted cuprates. We propose a resonant inelastic x-ray scattering 
experiment for Ni if -edge to confirm hole binding around the Ni impurity. 

PACS numbers: 74.72.-h, 71.10.Fd, 71.55.-1 



Atomic substitutions for copper in high-temperature 
cuprate superconductors induce significant impacts 
on macroscopic and local physical properties of the 
cuprates [ij. Divalent transition metal ions such as zinc 
and nickel have been frequently substituted for Cu. It is 
well known that nonmagnetic Zn impurities suppress the 
superconducting transition temperature more strongly 
than magnetic Ni impurities. Concerning magnetic prop- 
erties, neutron scattering experiments have shown that 
Zn impurities enhance the antiferromagnetic (AF) cor- 
relation length only slightly, whereas Ni impurities sup- 
press incommensurate peaks strongly and stabilize the 
Neel ordering [3| . Such contrasting behaviors have been 
discussed in connection with the difference of spin state 
of Zn2+ with spin S* = and Ni2+ with 5 = 1 For 
Ni substitution, however, it has been argued from several 
experiments 0, 0, d, d, that a Ni impurity attracts a 
hole, giving Ni'^+ with S = ^. More recently, it has been 
suggested from x-ray-absorption-fine-structure (XAFS) 
measurements that the Ni impurity forms Ni^+L state 
with 5* = i (L represents a ligand hole), i.e., Zhang-Rice 
(ZR) doublet state, binding a hole on neighboring oxygen 
orbitals [§]. The formation of the doublet may explain 
several experimental facts for Ni-substituted cuprates as 
discussed in Ref . Q : small moment of Ni impurity 0, [l^ , 
weak suppression of coherence peaks reduction of 
magnetic resonance energy [l2l | , enhancement of pseudo- 
gap energy [11], and so on. In spite of accumulating ex- 
perimental evidences of hole binding around Ni impurity, 
there is no theoretical support clarifying the formation of 
the ZR doublet as far as we know. 

In this Letter, we show that the ZR doublet is cer- 
tainly formed at Ni site embedded in the Cu02 plane, 
based on exact diagonalization calculations for small clus- 
ters with realistic parameter values. This theoretical re- 
sult together with experimental ones mentioned above 
puts restrictions on theoretical models of Ni-substituted 
cuprates. We propose a site-selective experimental tech- 
nique to confirm the presence of the ZR doublet, which is 
resonant inelastic x-ray scattering (RIXS) experiment for 
Ni K-edge. We predict low-energy structures inside the 



Mott gap, which can be direct evidence of hole binding 
around Ni impurity. 

For Cu and O ions, we consider Cu3(i2.2_j^2 and 02p„ 
orbitals. An additional orbital Sd^z^-r^ is included on 
Ni ion substituted for Cu. Apical oxygen 2pz orbitals 
below/above the Ni ion are also taken into account. By 
including hoppings between 3c? and 2p orbitals and the 
Coulomb interactions among 3d orbitals, the Hamilto- 
nian in the hole representation is given by Hdp = Hx+Hd 
with 
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where the summation of i runs over all of Cu and Ni 
sites, io denotes Ni sites, x (y) is the vector connecting 
neighboring Cu ions along the x {y) direction, io ± z rep- 
resents apical oxygen sites above/below Ni sites, and 6 
runs over x, y, and 2z. The operator dia {d[^^) is the 
annihilation operator for a 3d.j.2_y2 (3c?3^2„r2) hole with 
spin a at site i (io), Pi±i^ is the annihilation operator for 



a 2p hole at site i ± | with spin a 
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Tpd, Tp^, and T^'^ are the hopping 



integrals between "idx^-y^ and 2po-, between 3^322 .^2 and 
2pa, and between 3^322 _^2 and 2pz, respectively, with 

Tpd = VST^rf = ^T^d- The tetragonality of NiOe octa- 
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hedron is represented by introducing a, where < a < 1 
for elongated octahedron along the z direction. The pa- 
rameters Ed (ffp) and Ud are the energy level of Cu3dx2_y2 
{02p) and the Coulomb repulsion of Cu3d2,2_j^2 orbital, 
respectively, em, Um, U^-^, and i^Ni are the energy levels 
of Ni3(i orbitals, intra- and inter-orbital Coulomb repul- 
sions, and the exchange interaction, respectively, with 
U^i = C/^i -I- 2XNi- Hereafter we set Ed = 0. 

Following a procedure by Zhang and Rice [14], let us 
introduce three kinds of 02p Wannier orbitals, i.e., sym- 
metric, antisymmetric, and nonbonding. Among them, 
we can neglect the nonbonding orbital since it is decou- 
pled from other orbitals. For a system with N unit cells, 
the symmetric and antisymmetric Wannier orbitals are 
given by (t>f^ = -iN-2 Y,^^e'^'^Psi^iSxPx]i<7 ~ SyPyi^a), 
and^?. = iN'^ J2ke'^'Pak[V2SxSy{SyPxk„+SxPyi,a)- 
ia{Sl -f Sl)pzkcr], where S^(^y) = sin^^, /3sk = (5^ -I- 
S'y)-3, /3ak = (3sk{2SlSl + a^p-^)-^, and the operators 
PxkdT, PykCT, and Pzka are the Fourier transformations of 
Pi+^<y,Pi+ia, and ^(pi+zCT-Pi-ztr), respectively Using 
these Wannier orbitals, (1) reads 

+2V2^;<iE<v;i>j^+H-c. 

iojcr 
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where nj = ^ Ek /^Jk e''"^'"^^ ^(j = ^Ek/?.k(5^ - 
^2)e^k(i-j)^ ^ ^Y.kPlkPak^'^'^'-'^- Their 

value is the largest if i — j is small: for example, Ti^j = 
0.958, r/^j = 0, T/f^j+^ = -0.258, and r/^j = 0.915. In 
the case of a single Ni impurity in the system, it is con- 
venient to introduce a new antisymmetric operator given 
by = P-' Ej r[^i<t>!a ^ith = since 0?, 

only couples to the impurity. Then, the number of the 
orbitals in the system becomes 2N + 2. 

We use the Lanczos-type exact diagonalization tech- 
nique on a TV = y/S x -s/S-unit-cell cluster with a single Ni 
impurity under periodic boundary conditions [l5l |. The 
parameter values have been estimated to be T^dS^ 0.95 ~ 
1.3 eV, [/d ~ 8 - 10.5 eV, £p ~ 2 - 3.5 eV In the 
present study, we take T^d = 1 eV, Ud = C^Ni = 8 eV, 
i^Ni — 0.8 eV, and Ep = 3 eV. is an unknown param- 
eter and a is taken to be 1/ V2 for simplicity. 

Insulating cuprates are known to be charge-transfer 
(CT) type, where the CT energy A between Cu and O 
is smaller than Ud [17]. Nickel oxides also belong to the 
CT typejbut the charge transfer energy Ani is larger 
than A [13. A is given by A = E{d^^L) - E{(P) = 
Ep, where E is the energy of a Cu-0 unit for a given 
configuration without hopping terms. Similarly Ani — 
E{(fL) - E{d^) ^ A-Em- Um + ^K^- We assume that 
An! for the Ni impurity is similar to the bulk systems, 
since the CT energy is predominantly determined by the 
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FIG. 1: (a) Averaged hole number (n) on each orbital in 
undoped ground state of an 8-unit-cell cluster as a function 
of the CT energy difference A = Ani — A. Filled upward and 
downward triangles represent Ni 3d^2_y2 (denoted as d) and 
3^322 _,,2 (d'), respectively, (p" and on Ni site are shown 
by open triangles. Filled (open) circles denote Cu d (Cu (p") 
orbitals at the nearest neighbor (n.n.) sites from Ni. (b) The 
difference of (n) between one-hole-doped and undoped ground 
states. Open squares and diamonds denote (p" orbitals at 
the second and third n.n. sites, respectively. Other symbols 
are the same as (a). The vertical broken line represents the 
position where the ground-state symmetry changes. 



energy level of 3d orbitals when ligand and environment 
are the same Therefore, we take Ani > A in the 

present study of Ni impurity, leading to ejsh < — 5.6 eV. 

Figure [Ha) shows hole number (n) on each orbital of 
the cluster as a function of the difference of the CT ener- 
gies A = Ani — a in the undoped case (9 holes in the 8- 
unit-cell system) . We find that the total number of holes 
at the Ni site (the sum of Ni d, Ni d', Ni (j)'' and (p") is 
almost two (1.99 at A = 0.4 eV and 2.04 at A = 3.9 eV), 
which means that nominally Ni^+ is formed. The local 
spin state is of high spin {S — 1). Remaining holes are 
distributed on Cu d, Cu 0^ orbitals, forming Cu^"*" on 
each Cu site. Note that (n) on the Cu orbitals at the 
second and third nearest-neighbor (n.n.) sites from Ni is 
nearly the same as those of the n.n. sites. 

Next we examine how a hole introduced into the sys- 
tem is distributed among the orbitals. In the range of 
A shown in Fig. [U we find that with increasing A the 
ground state with one additional hole (10 holes in total) 
changes quantum number at A 1.2 eV from the total 
spin 5=1 and A2 irreducible representation of 6*4^ (de- 
generated with Bi) to S" = and Bi. The difference of 
hole number, S{n), on each orbital between the one-hole 
doped and undoped ground states is plotted in Fig.ljb). 
We find that for A > 1.2 a hole introduced in the system 
occupies predominantly on Ni 0" orbital, i.e., an oxygen 
Wannier orbital that mainly couples to Ni3c?a;2_y2. In 
total, more than 60% of the hole enters into Ni-related 
orbitals. This is consistent with experimental sugges- 
tions [1, 0, S S S IB that a Ni impurity may bind a hole 
and form the ZR doublet. Note that negative S{n) for 
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Ni d! is the consequence of strong energy gain by accom- 
modating holes into the orbitals with — symmetry. 
For A > 1.2, the doped hole is distributed mainly on 
oxygen orbitals away from the Ni impurity. Since an ef- 
fective CT energy for Ni oxides has been reported to be 
~ 5 eV it is reasonable to consider that A > 1.2 eV. 
Therefore, the present result strongly supports the bind- 
ing of a hole around a Ni impurity in real Ni-substituted 
cuprates. 

In order to confirm the stability of the ZR doublet, we 
compared binding energy of the ZR doublet on a NiOg 
cluster with that of the ZR singlet on a CUO4 cluster. 
We found that the ZR doublet is more stable than the 
ZR singlet when Ani is larger than A by 0.3 eV. This is 
qualitatively consistent with the data shown in Fig.[ljb), 
implying that the stability of the ZR doublet is governed 
by local character around Ni. In the limit of (C/^, C/ni, 
A, Ani) 3> Tpd that is approximately satisfied in the 
present parameter set, the binding-energy difference of 
the doublet and the singlet, (5i?B, may be given by [l§| 
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We notice that dominant negative contribution to SEb 
comes from the third term in (j4|). This implies the im- 
portance of the condition that J/ni (~8 eV) is not much 
larger than Ani (~5 eV). 

When a hole is bound to a Ni impurity, localized Cu 
spins are expected to be unaffected and thus AF correla- 
tion may remain. Defining 5*? = (nf^ — nf ^ + </>; | — 
l)/2, we show in Fig. [2] the equal-time spin-spin 
correlation function {S?S?) between two Cu sites in un- 

doped and one-hole doped cases with A = 1.4 and 0.9 
eV. Here, A is chosen near the boundary of two regions 
shown in Fig. [T]Jb). In the undoped case, the correlation 
on the same (different) sub-lattice is positive (negative), 
implying the presence of AF order. For A = 1.4 eV, 
{S?S?) in the one-hole doped case is very similar to the 
undoped case. This is a natural consequence of hole bind- 
ing around Ni impurity. On the other hand, the correla- 
tion is dramatically changed from the undoped case when 
A = 0.9 eV. This is due to the destruction of AF order 
caused by carrier motion. 

To directly observe hole binding around Ni impu- 
rity, we need to use site-selective probes. One of 
them would be XAFS that has been performed re- 
cently Q. As another probe, we propose a RIXS ex- 
periment for Ni X-edge, which can detect charge ex- 
citations related to the bound hole. In the Ni K- 
edge RIXS, the emission of a photon with a dipole 
transition between Ni 4p and Ni Is states occurs res- 
onantly by tuning the energy of incoming photon to 
Ni K absorption edge. In the intermediate state, we 
introduce Coulomb interaction between 3d and Is-core 
holes ;201, given by iJis-ad = Eio,7,<T,a' '^L^W" 
where n?^^ is the number operator of Is-core hole. 
By assuming that the 4p photo-electron enters into 




1 2 3 4 5 
pair number 



FIG. 2: (a) Schematic diagram for an 8-unit-cell cluster (the 
inside of dotted lines) where the filled and open circles de- 
note Ni and Cu sites, (b) Spin-spin correlation function for 
the pairs numbered in (a). The upward (downward) trian- 
gles show data for A = 1.4 (0.9) eV, and the open and filled 
symbols denote the undoped and one-hole-doped cases, re- 
spectively. 



the bottom of the Ap bands [20J, the RIXS spectrum 
for Ni K-edge as well as Cu i^-edge is expressed as 



where Z3k = Yl'i a 'pIp io-^L + -^.c. with the creation 
operator (plpia) Is-core hole (4p electron). ^| 
denotes summation over Ni (Cu) sites for Ni (Cu) K- 
edge. Kj(^-) is the wave vector of the incoming (out- 
going) photon with energy i^i(f), and Auj = uji — LOf. 
G~^{u!i) = Ui + iV - Hdp - His-3d - His^ip, where i?is,4p 
is composed of the energy separation £^^^"4^ between the 
Ni (Cu) Is level and the bottom of the 4p band, and F is 
the inverse of relaxation time in the intermediate state. 
|0) is the ground state with energy Eq, \f) is the final 
state of RIXS with energy Ef. We use F = 1 eV and 
Uc — 4: eV JJi]. I (Auj) is calculated by using a modified 
version of the conjugate-gradient method together with 
the Lanczos technique. 

Figure [3] shows Cu if -edge (upper panel) and Ni K- 
edge (lower panel) RIXS spectra in undoped and one- 
hole-doped cases for A = 1.4 eV. The spectra at zero 
momentum transfer (K f — = 0) are shown for Cu 
K-edge. In the undoped case, the edge of the CT gap 
is located at ~ 2.1 eV and ~ 3.2 eV for Cu and Ni, 
respectively. The difference comes from different CT en- 
ergies. This is qualitatively consistent with experimental 
observations for La2Cu04 (~ 2.2 eV) and La2Ni04 4 
eV) [2l|. Furthermore, recent RIXS experiments for Ni- 
substituted La2Cu04 also show similar behaviors [22*1. In 
Cu K-edge RIXS, the spectrum around 4.5-6 eV is asso- 
ciated with the excitations from bonding to antibonding 
states in CUO4 plaquette [1^, . A similar structure ap- 
pears in the energy region from 6 to 7 eV for Ni K-edge 
RIXS. 

The Ni K-edge spectrum is strongly affected by hole- 
doping as shown in the lower panel of Fig. [3l In par- 
ticular, new spectral structures appear at Aui ~ 1.35 eV 
and ^ 2.13 eV within the CT gap. Examining eigenstates 
generating the new structures, we find that the hole num- 
ber of Cu is larger than that in the ground state, while 
the hole number at the Ni site is smaller. This means that 
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FIG. 3: RIXS spectra for Cu (upper panel) and Ni (lower 
panel) if -edges in an 8-unit-cell cluster with a Ni impurity. 
The solid and broken lines denote the undoped and one-hole- 
doped cases, respectively. A = 1.4 eV. The 5-functions (the 
vertical lines) are convoluted with a Lorentzian broadening of 
0.2 eV. 

the ZR singlet on Cu are dominating in the correspond- 
ing excited states. Since the hole number at the Ni site is 
large in the ground state due to the formation of the ZR 
doublet, the new spectral structures are associated with 
excitations from the ZR doublet to ZR singlet. There- 
fore, if we observe spectral weights inside the CT gap for 
Ni-substituted cuprates, it can be identified as direct ev- 
idence of hole bounding around Ni impurity. Note that 
the spectrum for Cu X-edge is almost unchanged upon 
hole-doping although a small hump is seen around 1 eV. 

We have neglected direct 0-0 hoppings in ([1]). The 
hoppings may cause the decrease of A, leading to less 
stability of the ZR doublet as seen from However, 
we found that the boundary at A ~ 0.3 eV for single clus- 
ters without the hopping parameter Tpp increases only by 
0.1 eV for a reahstic value of Tpp {—0.5 eV). 



The formation of the ZR doublet at Ni site imposes re- 
strictions on modeling impurity effects in cuprate super- 
conductors, in particular, on the construction of single- 
band models like i-J-type model [H [H, [13, [H, [H, [sO] ■ 
We should treat a Ni impurity as a 5' = ^ spin in hole- 
doped cases. The spin is expected to couple to the neigh- 
boring Cu spins antiferromagnetically with an interaction 
parameter determined by competition between superex- 
change processes and the process of the motion of the (j)" 
hole. Other holes that hop to the impurity position may 
feel repulsive interaction. This model is close the those 
used in Ref. More precise effective i- J- type model 

is now under construction. 

In summary, we have carried out the numerically ex- 
act diagonalization on the eight-unit-cell cluster with a Ni 
impurity site representing Ni-substituted cuprate super- 
conductors. Using realistic parameters of the cuprates, 
we have found that a hole can be bound in the Ni04 pla- 
quettes forming the ZR doublet. This is due to the fact 
that the ZR doublet is more stable in energy than the ZR 
singlet on Cu. This finding supports theoretically recent 
experimental suggestion of the presence of the ZR dou- 
blet. Also this imposes strong restrictions on modeling 
Ni-substituted cuprates. We have proposed that the hole 
binding can be seen in RIXS for Ni K-edge. We hope 
that the proposed RIXS experiment will be done in the 
near future. 
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